High-Field Electrical Transport in Single- Wall Carbon Nanotubes 
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Using low-resistance electrical contacts, we have measured the intrinsic high-field transport prop- 
erties of metallic single-wall carbon nanotubes. Individual nanotubes appear to be able to carry 
currents with a density exceeding 10 9 A/cm 2 . As the bias voltage is increased, the conductance 
drops dramatically due to scattering of electrons. We show that the current-voltage characteris- 
tics can be explained by considering optical or zone-boundary phonon emission as the dominant 
scattering mechanism at high field. 

PACS numbers: 73.50.-h, 73.61.Wp, 73.50.Fq, 72.10.Di 



The potential electronic application of single-wall car- 
bon nanotubes (SWNTs) requires a detailed understand- 
ing of their fundamental electronic properties, which are 
particularly intriguing due to their one-dimensional (ID) 
nature (l). Metallic SWNTs have two ID subbands cross- 
ing at the Fermi energy. In the ideal case the resistance is 
thus predicted to be h/Ae 2 or 6.5 kf2. In early electrical 
transport experiments, however, the nanotubes typically 
formed a tunnel barrier of high resistance of ~1 Mf2 with 
the metal contacts ||||. Consequently the bias voltage 
dropped almost entirely across the contacts, and tunnel- 
ing dominated the transport. A number of interesting 
phenomena have been observed in this regime. At low 
temperatures, Coulomb blockade effects prevail ^M- At 
relatively high temperatures, the transport characteris- 
tics appear to be described by tunneling into the so-called 
Luttinger liquid — a unique correlated electronic state in 
ID conductors which is due to electron interactions 

One of the most important questions that remains to 
be addressed is how the electrons traverse the nanotubes, 
i.e., whether ballistically or being scattered by impurities 
or phonons. The unusual band structure of metallic tubes 
suggests a suppression of elastic backscattering of elec- 
trons by long-range disorder p| . Long mean free paths for 
electrons near the Fermi energy have indeed been inferred 
from regular Coulomb oscillations and coherent tunnel- 
ing at low temperatures dj7j]. However, there has been 
no transport study of electrons with significant excess en- 
ergy above the Fermi energy. It is not clear whether such 
electrons would experience strong scattering and what 
type of scattering mechanism would dominate. 

In this Letter we present electrical transport measure- 
ments of individual nanotubes using low-resistance con- 
tacts (LRCs). In contrast to the high- resistance contacts 
(HRCs), a bias voltage applied between two LRCs estab- 
lishes an electric field across the nanotube which accel- 
erates the electrons, enabling transport studies of high- 
energy electrons. We find that individual SWNTs can 
sustain a remarkably high current density of more than 
10 9 A/cm 2 . The current seems to saturate at high elec- 



tric field. We discuss possible scattering mechanisms and 
suggest that optical or zone-boundary phonon emission 
by high-energy electrons can explain the observed behav- 
ior. An analytic theory based on the Boltzmann equation 
is developed which includes both elastic scattering and 
phonon emission. The numerical calculations reproduce 
the experimental results remarkably well. 
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(b) High-resistance contacts 
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FIG. 1. Typical / and dl/dV vs V obtained using (a) low- 
and (b) high-resistance contacts. The inset to (a) shows a 
0.86 x 0.86 /im 2 AFM height image of a typical LRC sample, 
where the z range is 10 nm. The inset to (b) plots dl/dV vs 
V on a double-log scale for the HRC sample. 

The inset to Fig. 1(a) shows an atomic force micro- 
scope (AFM) image of our typical LRC sample. The 
20 nm thick, 250 nm wide Ti/Au electrodes are embedded 
in thermally-grown Si02 with a height difference of less 
than 1 nm which minimizes the deformation of the nan- 
otubes near the electrodes. This is achieved by electron- 
beam lithography and anisotropic reactive-ion etching of 
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SiC>2 using a single layer of PMMA as both electron-beam 
resist and etching mask, followed by metal evaporation 
and lift-off. The electrodes are cleaned thoroughly in 
fuming nitric acid. The nanotubes are then deposited 
on top of the electrodes from a suspension of SWNTs 
ultrosonically dispersed in dichloroethane. We find that 
brief annealing of the electrodes at 180°C improves the 
reproducibility of the contact resistance. Only nanotubes 
with apparent height of ~1 nm under AFM are chosen 
for transport measurements, which are presumably indi- 
vidual SWNTs. Metallic nanotubes are selected based on 
the absence of gate effect on transport at high temper- 
atures |^]. This procedure yields a typical two- terminal 
resistance of individual metallic tubes of less than 100 kSl 
(the lowest is 17 kf}) at room temperature, as compared 
to the ~ Mfi resistance using Pt as the contact material 
in previous experiments S . Similar reduction in contact 
resistance has also been achieved in a different contact ge- 
ometry ||To|| . The exact mechanism for the low contact re- 
sistance is unclear. However, clean flat single-crystalline 
gold facets may increase the coupling by increasing the 
effective contact length over which a small tube-electrode 
separation is realized [jll| . 

Figure 1 shows the typical two-terminal current / and 
differential conductance dl/dV vs voltage V obtained us- 
ing LRC (Au) and HRC (Pt). dl/dV is acquired simulta- 
neously using a standard ac lock-in technique. The room- 
temperature zero-bias resistance of the two samples are 
40 kil and 670 kf2 respectively. For both samples, the 
zero-bias conductance G decreases monotonically as the 
temperature T decreases Q. The large-bias- voltage de- 
pendence of dl/dV, however, is notably different. For 
the LRC sample, dl/dV increases with increasing bias, 
reaching a maximum at ^100 mV. As the bias increases 
further, dl/dV drops dramatically. In contrast, the HRC 
sample exhibits a monotonic increase of dl/dV as a func- 
tion of voltage up to 1 V. The inset to Fig. 1(b) plots 
dl/dV vs V on a double-logarithmic scale for the HRC 
sample, in which it appears that dl/dV can be fit with 
a power-law function for large bias. Both the tempera- 
ture dependence of G and the bias-voltage dependence 
of dl/dV for the HRC sample are typical of individual 
SWNTs and ropes with similar or lower conductance val- 
ues Q , which are attributed to the suppressed tunneling 
density of states in a Luttinger liquid fl^| . The similar 
behavior around zero bias for the LRC sample suggests 
that it comes from the same origin. In the remaining of 
the paper, we focus on the large-bias behavior for the 
LRC samples. 

We have further extended the I-V measurements up 
to 5 V as shown in Fig. 2. Strikingly, the I-V curves 
at large bias measured at different temperatures between 
4 K and room temperature essentially overlap with each 
other. The current at 5 V exceeds 20 fiA, which corre- 
sponds to a current density of more than 10 9 A/cm 2 if 
a spatial extent of the 7r-electron orbital of ^3 A is used 



to estimate the current-carrying cross section. From the 
shape of the I-V curves, it is clear that the trend of 
decreasing conductance continues to high bias. Extrapo- 
lating the measured I-V curves to higher voltage would 
lead to a current saturation, i.e., a vanishing conduc- 
tance. Interestingly, the saturation current seems to be 
independent of the distance between the electrodes |13] . 




-5 -4 -3 -2 -1 1 2 3 4 5 
V(V) 

FIG. 2. Large-bias I-V characteristics at different temper- 
atures using low-resistance contacts for a sample with an elec- 
trode spacing of 1 /im. The inset plots R = V/I vs V. 

We find that the resistance, R = V/I, can be fit re- 
markably well by a simple linear function of V over al- 
most the entire range of applied voltage (see Fig. 2 inset): 

R = R + V/I , (1) 

where Rq and Jo are constants. Only at high voltages 
near 5 V some samples start showing slight deviation 
from this linear behavior. From the slope of the linear 
part of R-V we find Iq ~ 25 /zA which is approximately 
the same for all of the samples we have measured. 

At first sight the current saturation might be explained 
from the band structure. Current in metallic nanotubes 
is carried by two propagating ID subbands. In the ab- 
sence of scattering, the chemical potentials of the right 
and left moving states will differ by the applied voltage 
eV. At low voltages this leads to an Ohmic response, 
but when eV exceeds the Fermi energy of the ID sub- 
bands, the left moving states will be completely depleted 
and the current will saturate. The Fermi energy, mea- 
sured relative to the nearest band edge, is approximately 
2.9 eV. Experimentally, however, the current starts sat- 
urating at a much lower voltage. An alternative model is 
thus needed to explain the saturation. 

We expect the measured resistance to be a combination 
of the resistance due to the contacts and the resistance 
due to backscattering along the length of the nanotube. 
The current saturation is unlikely to arise from an in- 
creased contact resistance at high voltages since the con- 
tacts would then behave like high-resistance tunneling 
contacts, and one would expect to see features in the I- 
V associated with tunneling into the ID subbands of the 
nanotube. The measured I-V, however, is featureless. 
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We thus focus on the effect of backscattering in the 
nanotube. The behavior of R-V suggests that in ad- 
dition to a constant scattering term, which most proba- 
bly comes from contact scattering or impurity scattering, 
there is a dominant scattering mechanism with a mean 
free path (mfp) which scales inversely with the voltage. 

Electrons can backscatter off phonons and other elec- 
trons. Electron-electron scattering is appealing at first, 
since it does not involve heating the lattice. The only 
electron-electron scattering that contributes to resistiv- 
ity is Umklapp scattering |Q with a scattering rate di- 
rectly proportional to the electron temperature T c . This 
gives V/I ~ T e , T c will be determined by how fast the 
heat can escape from the tube. If we assume that all the 
heat produced is carried by electrons into the leads and 
that the temperature along the tube is uniform [jl5j, we 
have IV = 4(7r 2 /3)(fcBT' c ) 2 //i, where the left-hand side is 
the rate at which heat is produced, and the right-hand 
side is the heat current carried by the two ID channels. 
Hence we expect that I ~ V 1 ^ 3 , which cannot describe 
the experiments. We have verified this V 1 ! 3 behavior 
by numerically solving a Boltzmann equation similar to 
that discussed below. Luttinger-liquid effects, which have 
been ignored in the above arguments, tend to enhance 
Umklapp scattering at low energies Q. This would 
make the agreement with experiment even worse. 

This suggests that we must consider scattering from 
phonons. The coupling will be strongest for phonons 
which compress and stretch bonds on the lattice scale. 
There are three possible categories of phonons: (1) 
twistons or long-wavelength acoustic phonons jlq], (2) 
optical phonons which are derived from the in-plane E2 92 
mode of graphite with a frequency of 1580 cm -1 , and (3) 
in-plane zone-boundary phonons with momentum which 
connects the two Fermi points of graphene. While zone- 
boundary phonons are not directly observable optically, 
force-constant models put their frequency in the range 
1000 — 1500 cm -1 in graphene [0. Twiston scattering 
is most likely not relevant since the scattering rate is 
smaller than the optical or zone-boundary phonon scat- 
tering by ~ T/Qjj, where T is the lattice temperature 
and 0£> ~ 2000 K is the Debye temperature It is 

unlikely that the lattice temperature could be that high. 
Moreover, twistons may be pinned by the substrate. 

Now we discuss backscattering due to the emission of 
optical or zone-boundary phonons. A related effect has 
been discussed previously in the context of semiconduc- 
tors ]lq| . The key point is that for an electron with en- 
ergy E to emit a phonon of energy HQ, there must be an 
available state to scatter into at energy E — hil. In the 
presence of an electric field £, electrons are accelerated, 
hk = e£. It is simplest to consider the case in which the 
coupling to the phonon is so strong that, once an electron 
reaches the threshold for phonon emission, it is immedi- 
ately backscattered. As indicated in the schematic in the 
inset to Fig. 3, a steady state population is then estab- 



lished in which the right moving electrons are populated 
to an energy hft higher than the left moving ones. The 
current carried in this state can be computed from a Lan- 
dauer type argument to be 

J = {Ae/h)Ml. (2) 

If we choose Ml = 0.16 eV (corresponding to 1300 cm -1 ), 
this leads to a saturation current of 25 fiA, which is inde- 
pendent of sample length and agrees very well with the 
measured saturation current. 

In this picture, the mfp for backscattering phonons £q 
is equal to the distance an electron must travel to be ac- 
celerated to an energy above the phonon energy: £q = 
HQ/e£. This may be combined with a constant elastic 
scattering term via Mathiesson's rule to obtain an effec- 
tive mfp, £~g = ^<T 1 + ^ 1 , where £ c is the elastic scatter- 
ing mfp. The resulting resistance, R = (/i/4e 2 )(L/£ e ff), 
then has the empirically observed form of Eq. (1) with 
Ro = 0/4e 2 )L/4 and J given in Eq. (2). 

To put the above interpretation on a more quantitative 
basis, we consider the Boltzmann equation for the distri- 
bution functions fL,n{Ek,x,t) of left and right moving 
(L, R) electrons. Details will be provided in a future 
publication. 

[d t ± v F d x ± v F e£d E ] fh.R = [dt/i,ii] col ■ ( 3 ) 

Here vf is the Fermi velocity, and we have chosen to 
express the momentum dependence of / in terms of 
Ek = ±TiVFk. The left-hand side describes the col- 
lisionless evolution of the electrons in the presence of 
an electric field £. For the collision term on the right, 
we consider a sum of three terms: (1) Elastic scatter- 
ing, [d t fh\ c = {vf/QUr - h), where £ c is the elas- 
tic mfp. (2) Backscattering from phonons, [<3t/rJ p b = 

(v F /£ P b) [(1 - f L )ft - h(l - f R )] ■ Here /± are eval- 
uated at E ± HQ,. £ p b, which depends on the strength 
of the electron-phonon coupling, is the distance an elec- 
tron travels before backscattering once the phonon emis- 
sion threshold has been reached. This should be con- 
trasted with £q , the distance required to reach the thresh- 
old. We assume that the phonon temperature is much 
less than the phonon energy of ^2000 K, so that the 
Bose occupation factors can be ignored. Finally, we con- 
sider (3) forward scattering from phonons, [<9t/L] p f = 

(v F /£ p{ )[(i-fL)fi-f L (i-fE)]- 

The effects of the contact resistance may be included 
as a boundary condition at the ends of the tube. For 
instance at the left contact [x = 0) we have, 

f R (E, 0) = t\ f (E - ml) + (1 - tl)f L (E, 0), (4) 

where t 2 L is the transmission probability for the contact, 
and/ (£-/ii) = (exp[(S-/i L )/fc B J 1 ] + l)~ 1 is the Fermi 
distribution function of the left contact with electrochem- 
ical potential [II and temperature T. 
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FIG. 3. Numerical calculation of I-V characteristic by solv- 
ing Boltzmann tranport equation including elastic impurity 
scattering and phonon emission as depicted in the left inset. 
See text for the parameters used. 



that all of the heat could be transmitted through the 
contacts. However, the nanotube is in intimate contact 
along its entire length with the substrate, which may 
be regarded as a thermal reservoir at T < 300 K. It 
would clearly be desirable to study further the nature 
of the thermal contact between the nanotube and sub- 
strate. Measurements on suspended nanotubes may pro- 
vide some useful information. 

We thank R.E. Smalley and coworkers for providing 
the indispensable nanotube materials, M.P. Anantram, 
H. Postma and S.J. Tans for discussions, and A. van den 
Enden for technical assistance. The work at Delft was 
supported by the FOM and the work at Pennsylvania by 
the NSF under grant DMR 96-32598. 



We have solved Eqs. (3,4) numerically to obtain the 
steady state distribution function Jl,r{E, x) in the pres- 
ence of an applied voltage jj,L — /Ar = eV as a function of 
hfl, £ e , ^ p b,pf, and t\ R . The current is then simply given 
by I = (4e//i) / dE(fi, — /#). Figure 3 shows an exam- 
ple of the numerical calculation of I-V characteristic for 
a sample length of L = 1 /im. The parameters used in 
the plot are Titt = 0.15 eV, t\ R = 0.5, t e = 300 nm, 
£pb = 10 nm, and £ p { = oo. The resemblance to the 
experiment is remarkable. It is interesting to note that 
the current is insensitive to the contact scattering for 
V > 0.5 V. Contact scattering affects only the low bias 
resistance, giving rise to the positive curvature in R-V 
near 1^ = 0. 

Assuming local thermal equilibrium the Boltzmann 
equation may be used to derive hydrodynamic equations 
which govern the transport of charge and energy. These 
equations may then be solved analytically and give re- 
sults which agree well with the simulations. They show 
that: (i) The empirical formula [Eq. (1)] is exact in 
the limit eV <C MlL/£ p ^,, which means that the energy 
gained by an electron within distance £ p b must be much 
less than the phonon energy, or equivalently, i? p b <C in- 
(ii) For larger V the simple formula breaks down, and 
in the limit of very large V , the resistance becomes con- 
stant, R = (/i/4e 2 )i(^ ( T 1 +^pb)- For the parameters used 
in Fig. 3, the crossover voltage is roughly 15 eV. Indeed, 
there appears a small negative curvature at 5 V in the 
V/I vs V plot (inset Fig. 3), which signals the beginning 
of the breakdown of the empirical formula. The curva- 
ture would be less pronounced if a shorter value for £ p b 
is used. We note that 10 nm seems rather short. An es- 
timate using a simple deformation potential model gives 
£ p b ~ 150 nm for a nearly armchair nanotube. More 
work is needed to have a more accurate estimate of the 
electron-phonon coupling strength. 

We have assumed that the heat generated in the tube 
escapes sufficiently quickly to avoid raising the lattice 
temperature too high. A simple estimate of the nan- 
otube's thermal conductivity indicates that it is unlikely 
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